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DIAMAGNETISM VERSUS ROBIN CONDITION AND CONCENTRATION 

OF GROUND STATES 

AYMAN KACHMAR 


Abstract. We estimate the ground state energy for the magnetic Laplacian with a Robin 
condition. In a special asymptotic limit, we find that the magnetic field does not contribute to 
the two-term expansion of the ground state energy, thereby proving that the Robin condition 
weakens diamagnetism. We discuss a semi-classical version of the operator and prove that the 
ground states concentrate near the boundary points of maximal curvature. 

1. Introduction 

This paper is motivated by two different questions. The first question concerns the analysis of 
the ground state energy and the concentration of the ground states for the magnetic Laplacian 
with a Robin condition and a semi-classical parameter, and is a continuation of the work in 
ia uni. The second question is around the influence of the Robin condition on diamagnetism. 
We will find that these two questions are intimately related and we will get satisfactory answers 
for both. Besides the concentration of the ground states near the points of maximal curvature, we 
will identify the optimal strength of the Robin condition/magnetic held such that diamagnetism 
occurs to leading order of the energy. 

The results in this paper are valid in an open set 17 C R 2 . We will assume that the boundary 
r = <917 of 17 is C 3 smooth, compact and consists of a finite number of connected components. 
Our assumptions allow for 17 to be an interior or exterior domain, and the smoothness of the 
boundary ensures the existence of a normal vector every where on the boundary. We will denote 
by v the unit outward normal vector (held) of <917. 

1.1. Concentration of ground states. Let us introduce the semi-classical magnetic Laplacian 
that we will study. Consider the magnetic potential: 

M 2 9 {xi,x 2 ) eA A 0 (xi,x 2 ) = (-£2,0). (I l l) 

This magnetic potential generates the constant magnetic held: 

B : = curl Ao = 1. (1.1.2) 

We are interested in the same magnetic Laplacian studied in [10] which involves four parameters, 
the strength of the magnetic held b > 0, the semi-classical parameter h > 0, two parameters 
7 € M \ {0} and a € M that will serve in defining the boundary condition. The operator is 

CV; n = -(W - ibA 0 ) 2 in L 2 ( 17), (1.1.3) 

with a boundary condition of the third type (Robin condition) 

v • (/iV — iAo)u + h a "fu = 0 on <917 . (1.1.4) 

This operator can be defined via Friedrich’s Theorem and the closed semi-bounded quadratic 
form, 

Qh,b,M '■= II ( W “ ibA o)u\\l^n) + / i 1+ "7 [ \u(x)\ 2 dx. (1.1.5) 

■J dQ. 

This quadratic form is defined in the ‘magnetic’ Sobolev space 

^-W°) = {«€L 2 (17) : (V-irtAo)€L 2 (17)}. (1.1.6) 
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The parameters a and 7 serve in controlling the ‘strength’ of the boundary condition in (11.1.51) . 
As we shall see, the sign of 7 and the values of a have a strong influence on the spectrum of the 
magnetic Laplacian in (II. 1.31) . Notice that 7 = 0 corresponds to the extensively studied magnetic 
Laplacian with Neumann condition (cf. [5j|2]), while 6 = 0 corresponds to the Laplacian without 
a magnetic field. That justifies the assumption 7 / 0 and 6 > 0. 

The ground state energy (lowest eigenvalue) of the operator in (11.1.31) is: 


6 , a, 7 ) 


ieHl 


inf 


O "’ 7 

^h,b,Q 


(u) 


'k 0 ( n )\{o} \\ u \r L 2 m 


(1.1.7) 


We will study the asymptotic limit where the semi-classical parameter h tends to 0_|_, while the 
parameters 6 , a and 7 are assumed fixed. In this regime, we can reduce to the case 6=1 simply 
by observing that, for all 6 > 0 , 

Hi(h; 6 , a , 7 ) = 6 2 /Ui ( 6 - 1 h; 6 = 1 , a, 6 ~ 1 + " 7 ) , ( 1 . 1 . 8 ) 

and that as long as we assume 6 fixed, 6 _1 /i —>• 0 when h —>• 0 +. 

The results in m distinguish between two situations. The first one corresponds to a > 1/2 
and is fairly understood: A two term asymptotic expansion of the ground state energy in (11.1.71) 
is established; the ground state energy is in the discrete spectrum (cf. |S]); and the ground states 
are localized near the boundary points where the curvature is maximal. 

The second situation corresponds to a < ^ and is less understood. Here the sign of 7 will 
play a dominant role. The contribution in this paper will clarify the situation when 7 < 0. For 
7 > 0 , the ground state energy satisfies 


/ii(/i; 6 , a, 7 ) = bh + ho(l) (h —>• 0 +). 

For 7 < 0, the behavior of the ground state energy is completely different and displayed as 
follows, 

Hi (h; 6 , a, 7 ) = — 7 2 /i 2a + h 2a o( 1 ) (h -A- 0 +). 

Note that this asymptotic expansion does not involve the strength of the magnetic field 6 . Again, 
the ground state energy is an eigenvalue, as long as the semi-classical parameter h is sufficiently 
small (cf. [ 8 ]). In |10| . it is proved that the ground states concentrate near the boundary (when 
a < I; and 7 < 0). The natural question is then weather one can refine the concentration near 
some special boundary points, e.g. points of maximal curvature. We will give an affirmative 
answer to this question in Theorem 11.11 below. 

Since the boundary is assumed smooth, there exists a geometric constant 1 0 € (0,1) such 
that, if dist(x, <912) < to, then we may assign a unique point p(x) € dQ such that dist(x,p(x)) = 
dist(x, (912). The function k(-) denotes the curvature along the boundary. 

In the statement of Theorem 11.11 re max is the maximum of the curvature along the boundary, 

c = V- 2 “ • 

7 

n € N is the smallest positive integer satisfying 

. , 1 — 2 a 1 

(n + 1--> 7 , (1-1-9) 

1 — a 2 

and e„(C) is the auantitv that we will introduce in (12. 2. 51) below. As C ^ 0 _|_, e n (£) behaves like 
(cf- Remark 12.61) . 

Now we are ready to state: 

Theorem 1.1. Suppose that 6 > 0, a < ^ and 7 < 0. 
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(1) As h —> 0+, the ground state energy in (II. 1.71) satisfies 

- 7 2 h 2a + e n (Ch 2 h 2a + 7 % a x h 1+a + h 1+a o( 1 ) if | < a < \ , 


Ri(h]b,a, 7 ) = < 


if a = ^ , 


-7 2 /l 2/3 + (4 ^b 2 + 7 K ma x)/l 4/3 + /l 4/3 o(l) 

„ -7 2 h 2a + 7 «max/i 1+Q + h 1+a o( 1 ) 

(2) There exist constants p G (0, ^), if G (0, -j), C > 0 and /ig > 0 such that, for all 


if a < | 


h G (0,/iq), ewery ground state Uh of (II.1.71) satisfies 


/ 




' ^bnd 

where 


\uh£f dx < Cex p —/i 2 ( 1 -“) 


) [ K,d" 


°~2 


dx < Cexp —h 2< - 1 ~ a '> 


flint = {x G fl : dist(x, dfl) > h 2 ( 1 -“)} and 

- ? 7 * 

flbnd — Gfl \ flint ■ ^max ft(p(x)) — ^ 2( 1 } 


Notice that the asymptotic expansions for b, a, 7 ) are compatible in the cases a = | and 
a < Formally, we get the expansion for a < | by taking 7 —>• —00 in the case a = f. 
Theorem 11.11 adds two improvements to the results in |10| by 

• establishing a two-term expansion for the ground state energy; 

• refining the concentration of the ground states near the points of maximal curvature. 

Remark 1.2. The magnetic field is assumed constant in Theorem 11.11 but the methods in the 
this paper should allow for dealing with non-constant C 1 magnetic fields. 


1.2. Diamagnetism. Here we will discuss the question of diamagnetism. We will find that 
imposing a Robin condition may slow diamagnetism (and even neglect this effect). Let j3 G K, 
H > 0 and C 6 (H) be the self-adjoint operator in L 2 (fl), 

CP{H) = -(V - iHA 0 ) 2 (1.2.1) 

with domain 

D(£P(H)) ={«G Hh Ao (H) : (V - iHA 0 ) 2 G L 2 (Q) and u • (V - iHA 0 )u + fu = 0 on dQ} . 

Here the magnetic Sobolev space Lf^ Ao (fl) is introduced in (II. 1.61) . Note that, for H = 0, £&( 0) 
is the Robin Laplacian, while for H > 0, C^{H) is the magnetic Laplacian with a (magnetic) 
Robin condition. Let <r(£^(iL)) be the spectrum of the operator £ 8 (H). We introduce the 
ground state energy, 

ih(/3-,H) = mf a (£P(H)). (1.2.2) 

The diamagnetic inequality yields, for all (3 G R and H > 0, 

- fh(/3;0) > 0. (1.2.3) 

In physical terms, this inequality refers to diamagnetism. It simply says that introducing a 
magnetic field increases the ground state energy. We will see that, when (3 —>■ — 00 , diamag¬ 
netism is weak in the sense that the difference pi (f3;H) — /2i(/3;0) is small. This property is 
a unique feature for the Robin condition as it fails for the Neumann and Dirichelt conditions. 
On the contrary, in simply connected domains, a Neumann boundary condition induces strong 
diamagnetism (cf. |2|). 

The asymptotic analysis of the spectrum of the Robin Laplacian is studied in many papers, 
cf- [HE, 0 EEi 13, HJ. In particular, as (3 —>• — 00 , the ground state energy satisfies, 

£i(/3;0) = -/3 2 + /3K max + /3o(l ). (1-2-4) 

We will write an asymptotic expansion for the magnetic ground state energy fi(f3', H ) valid when 
f3 —>• —00 and H —> 00 simultaneously. This is the content of Theorem 11.31 below. I 11 particular, 
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we will get a fair knowledge about the difference in 111.2.31 which measures the strength of 
diamagnetism. 

The statement of Theorem 11.31 requires a real-valued function ©(■) introduced in |1()| . For 
7 G R, define the ground state energy 


( 


@( 7 ) = inf 
£eR 


inf 

ueB 1 ^ R+)\{o} 

V 


[ ( W(t )\ 2 + \{t - 0«| 2 ) dt - 7k(0)| 2 ^ 

Jo _ 

11 112 > 


(1.2.5) 


where 

B\R + ) = {u e l 2 (r + ) : u',(t-£)ueL 2 ( R+)}. 

In jl'O] . it is proved that ©(•) is smooth, increasing and @( 7 ) > — 'y 2 for all 7 < 0. 
Now we are ready to state: 


Theorem 1.3. Let a € R \ {1} ; 0 < ci < C 2 and /5o < 0. Suppose that 

/3 < (3o and < H < C 2 1/3| • 

(1) If a > the ground state energy satisfies, as fi —» — 00 , 

fii(fcH) = <d(0)H + Ho{l). 

where ©o € ( 0 , 1 ) is a universal constant. 

(2) If a = 77 , the ground state energy satisfies, as (3 -A — 00 , 

H ) = H@ g (PH- 1 ' 2 ) + Ho( 1). 

(3) If | < a < 2 , the ground state energy satisfies, as ft —» —oo, 

H) = — ft 2 + e n{H(3 2 )/? 2 + fiKma.x + fio( 1 ) , 

where n is the smallest positive integer satisfying (11.1.91 and e n (-) is introduced in (|2.2.5p . 

(4) If a = Tj, the ground state energy satisfies, as (3 -A — 00 , 

H) = ^^— _ + K masj + (3o( 1 ). 

(5) If a < the ground state energy satisfies, as (3 -A — 00 , 

l(/3; H) = ~/3 2 + «max/3 + fio{ 1 ). 

Theorem II. 31 suggests that, in the limit j3 —>• — 00 , diamagnetism occurs to leading order when 
the strength of the magnetic field satisfies H ~ (3 a and a > 2. 

In the situation where H ~ /3 a and a < 2, diamagnetism occurs as a correction term and 
will compete with the correction term coming from the curvature of the boundary. According to 
Theorem 11.31 

• If | < a < 2, then diamagnetism occurs in the second correction term while the influence 
of the curvature occurs in the third correction term; 

• If a = |, both diamagnetism and curvature corrections appear in the second correction 
term; 

• If cr < dia-magnetism is weak and its contribution is negligible compared to the 

contribution of the curvature correction term (compare with (|1.2.41 1. 

Through this paper, the following notation will be used. C denotes a constant independent 
from the semi-classical parameter h. 0{h °°) is a quantity satisfying that, for all N £ N, there 
exist two constants ho € (0,1) and Cat > 0 such that, for all h € (0, ho), \0(h°°)\ < C n h N . 

The paper is organized as follows. In Section^ we analyze three auxiliary differential operators 
useful to prove Theorem 11.11 The proof of Theorem 11.11 occupies all of Section [3J Finally, in 
Section [4] we explain how to get the result in Theorem 11.31 from the existing results on the 
semi-classical magnetic Laplacian with a Robin condition, in particular those in Theorem 11.11 














DIAMAGNETISM VS ROBIN CONDITION 


5 


2. Analysis of auxiliary operators 

2.1. ID Laplacian on the half line. Here we introduce a simple ID operator that will play a 
fundamental role in the next sections. This operator arises naturally in the analysis of the Robin 
Laplacian without magnetic held (cf. [T2l I?]). The operator is 

d 2 

"Ho,o : = — ^2 in -k 2 (®4-) (2.1.1) 

with domain 

(ttGF 2 (M + ) : «'(0)=-tt(0)}. (2.1.2) 

The spectrum of this operator is { — 1} U [0,oo), and —1 is a simple eigenvalue with the L 2 
normalized eigenfunction 

uq(t) = \P1 exp(— t) . (2.1.3) 

2.2. Harmonic oscillator on the half-line. The key element in the proof of Theorem 13. II is 
the analysis of the harmonic oscillator 

^harm(C,0 = -^ + (Ct-0 2 hi L 2 (R + ) , (2.2.1) 

with domain 

D(n harm (C,0) ={neH\R + ) ■ T k u€L 2 (R + ), ke{ 1,2}, u'(0) =-u(0)}. 

Here £ > 0 and £ £ M are two parameters. Let us denote by (A n (%harm(C> £))) the increasing 
sequence of the eigenvalues of %harm(C>£) counting multiplicities. We will study the aymptotic 
behavior of the eigenvalue 

Al(L4arm(C,0) = inf (f 7 ^ harm(C)O)^ > (2.2.2) 

cis —y 0. 

By comparison with the operator in (I2.1.ip . we get: 

Lemma 2.1. For all £ > 0 and £ € M, it holds, 

Ai(74 arm(C? 0) > and A 2 (^harm(C ) 0) > 0 • 

The lower bound in Lemma IQ can be improved as follows: 

Lemma 2.2. There exists a universal constant Aq > 0 such that, if 0 < £ < 1 and |£| > AqQ, 
then 

Al(^harm(C,0)>-l + ^C 2 - 

Proof. Let u be an L 2 normalized ground state of the operator %harm(Ci£)- Let us write, 

Al(?4arm(C>0) = j (V| 2 + KC*"- 0 U \ 2 ) dr ~ l«(°)| 2 

( /•oo \ / roo 

J \u'\ 2 dr - |u(0)| 2 J +C 2 (y {W\ 2 + \(T-C 1 0 u \ 2dr ~\ u (°)\ 2 

We know that the lowest eigenvalue of the operator in (12.1.11) is —1. Let fi(A) be the eigenvalue 
of the operator 

--^ 2 +(t-A) 2 in L 2 (M + ), 

with the boundary condition tt^O) = — it(0). 

Now, it results from the min-max principle that 

Ai(^ harm (C,£)) > -(1 - C 2 ) + C 2 m(C _1 0 ■ 
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When |£ 1 ^| is sufficiently large, we get m(C 1 ^) > which is a result of the following two facts 
proved in [ 10 ]. 

lim fi(A) = oo and lim n(A) = 1 . 

A^—oo A—>-oo 


□ 


We will prove that: 


Theorem 2.3. Let n € N. There exist C > 0, a collection of vectors 

{m .7 = (Mj, 1 ) Mj, 2 ? ) Mj, 2 j+l) € M jy ^ , 

and a collection of vector functions, 

{Uj = (u jt i,u jj2 ,--- ,Uj, 2 j+l) € («S(K+)) 2j+1 j 

such that, if 0 < f < 1 and |£| < 1, then 

||(«hann(C,0 ~ A„K|| l2(]k+) < C^ 

where 


W n = 


n 2 j 

“o+EE^" p « p “i.p+>. 

j=i p=o 

uq is the eigenfunction in (I2.1.3p . and 


n 2j 

j'=i p=o 

Furthermore, 

hi = (mi,i = 2^1,2 = — 1 , A*i ,3 = l) • 

Proof. 

Step 1: Construction of (mi,ui). 

Here we construct Ml = (mi,1j Mi, 2 )Mi, 3 ) £ ® 1 3 and Hi = (111,1,111,2,111,3) such that the conclu¬ 
sion of Theorem 12.31 is valid for n = 1. 

Let us define 


Ai = —1 + Mi,iC 2 + Mi, 2 C£ + Mi,3^ 2 and wi(r) = u 0 (t) + CV,i(t) + C£hi,2(t) + £ 2 ui i3 (t). 

For simplicity of the notation, we will write H = 'Hharm(C, £)• Notice that, since + 1 ^ uq = 

0 , 


{H — Ai)u>i = ( 2 



Ul,l + (t 2 - 


+ e 2 


hl,l)u 0 


+ C£ 



,(yl^ + 1 ) ”'- 2_ 

111,3 + (1 - Mi,3)110 


(2 r + Mi,2)110 
+ f?i. ( 2 . 2 . 3 ) 


The remainder f?i is 

Ri =C 4 (t 2 - Mi,i)ni,i + C 3 £[(t 2 - Mi,i)ni,2 + (Ml,2 - 2 r)rxi >2 ] 
+ C 2 ^ 2 [(1 — Ml,l)ni ,3 + (Ml,2 — 2t )u \ t 2 + (1 — Ml,3)111,1] 
+ C £ 3 [(1 — Mi,3)111,2 + (mi,2 - 2 rK s] + £ 4 (1 — Mi,3)111 ,3 


(2.2.4) 
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We choose the coefficients and the functions in (12.2.311 so that all the terms on the left side vanish. 
This is possible since the operator —jjrpi +1 can be inverted in the orthogonal complement of the 
eigenfunction uq- That way we choose, 


f°° 1 f d 2 . 

Mi.i = J t 2 |u 0 (t )| 2 dr = -, = + {( r 2 _ Ml l ) Uo } 

r°° f d 2 \ -1 

Mi ,2 = - J 2r|tt 0 (r)| 2 dr = -1, ^ 2 + lj {(2r + Mi,i)^o} 

f°° ( d 2 \ _1 

Mi,3 = J \uo(r)\ 2 dr = 1, w lj3 = -(-—+ 1J {(1 - Mi, 3 )“o} • 


The operator + 1^ respects the Schwartz space <S(R+). The proof of this is standard 

(cf. |3], Lemma A.5]). Now, since |£| < AC, we infer from (12.2.3 1) and (12.2.4j) . 

\\(H — A 1 )u> 1 ||£2 (r+) < C(C 4 + £ 4 ) ■ 

Step 2: The iteration process. 

Suppose that we have constructed (Mj)j=i and (iij)” =1 such that 


(H Xn^Wn — R n T fn,(,£ • 


R n has the form 


Rr 


2n+2 

E c 2 n+ 2 ~ p evn, P 


p =0 


for a collection (u np ) of Schwartz functions that do not depend on ( and £, and the function 
/ n £ £ satisfies, 

ll/n,C,dlL 2 (R + ) < C (( 2n+4 + £~ n+i ) , 
where C is a constant independent of £ and £. 

We outline the construction of 


Mn+l — (Mra+1,1; Mra+1,2? ' j Mn+l,2n+3) € 


p 2?7/—|—3 


and ttn+l — (lln+1,1 j Hn+1,2> ‘ ' ' > Hn+l,2n+3) 


such that 


has the form 


(H A n+ i)u; n _(_i — Rn +1 T /n+i,C,£ j 

2n+4 

i? n+1 = EC 2 n+ 4 ’ P ^n + l, P , 

P=0 


for a collection (u n _|_i iP ) of Schwartz functions that do not depend on ( and £, and / n +i,c,£ 
satisfies, 

IIWdlL 2 (M + )<C(C 2 ” +6 + C 2 ” +6 ), 

where C is a constant independent of ( and £. 

We expand (H — A n +i)u; n -|_i and rearrange the terms in the form, 


2n+2 

(H A n _|_i)'?c n 4 -i = ^ ' C 

p=0 


2n+2- 


P£P 



^n+l,p+l T Vn,p 


Mn+l,p+l n 0 


2n+4 2n+6 

+ E C 2 n+ 4 - P C^n + l, P + E C 2 n+ 6 - P £ P Mn+l, P , 

p=0 p=0 

where the functions u n +i jP and M n +i,p are expressed in terms of the functions Uj tQ and the real 
numbers Hj, q . 
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All what we have to do now is to select the functions n n +i )P +i and the real numbers /i n+ i )P +i 
such that 

2n+2 ( / d 2 \ 1 

C 2n+2 P i P \ ( + -*- ) u n+l,p+l + Vn,p — Rn+l,p+lU0 ? = 0 . 

p=0 l V T / J 

To that end, we select /r n +i iP +i such that, 

POO 

Abi+l,p+l — / 'Un,pV'0dT , 

Jo 

so that 

Vn,p d"n+ l,p+1^0 -L hr L (1K_)_) . 

Finally, we define the function u n . i-i,p+i as follows, 


^n+l,p+l — ( df 2 ^ ) (fi n tP Mre+l,p+l^o) • 


□ 


As a consequence of Theorem 12.31 Lemma l2.1l and the spectral theorem, we get: 

Corollary 2.4. Letn £ N and A > 0. If |£| < AQ, then as £ —>• 0+, fde eigenvalue Ai("Hharm(C) 0) 
satisfies, 

n 2j 

Al(«hann(C,0) = -1 + E E /La-iC 2 -' ^ + 0(C 2n+2 ). 

i=i p=o 

Definition 2.5. Let n £ N and Ao &e fde universal constant in Lemma \2.2l For all £ (0,1), 
we define the following quantity 


e n (C) = inf{/n(C,£) : ICI < Cmax(A 0 ,1)} , 


where 


(2.2.5) 

( 2 . 2 . 6 ) 


72 2?2 

J=2p=0 

and (dj, P +1 ) are constants in Theorem \ 2. 31 

Remark 2 . 6 . Note that for n = 1 , /i(C> £) = ^C 2 — C£ + £ 2 j and 

min/i(C,0 = A(C,|C) = ^C 2 - 

Consequently, for all n £ N, as C —> 0+, 

en(C) = ^C 2 + 0 (C 4 )- 

2.3. A family of operators in a weighted space. Here we will study an operator that arises 
in many papers concerned with the semi-classical magnetic Laplacian (cf. [5l{TO|). Let h £ (0,1), 
C € (0,1), £ € R, fi £ M, 5 £ (0,1/2), m > 0, a £ (0,1), M > 0 and \fi\h 2~ 5 < 


Consider the Hilbert-space 

L 2 ((0,h- s ),adT ), a = l-(J3 + mh°)h 1 ' 2 T, INI L ^ h s Vadr) 


' rh~S \ >/ 2 

L M2 H 
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and the self-adjoint operator 

,M,h = ~ a _ 1 <9 r (a<9 r ) + (1 + h 1/2 A PtT ) (Cr(l - ^/3h 1/2 r) - 

= --^ + ((T-0 2 + (P + rnh")h 1 / 2 (l-(P + mh‘ T )h 1 / 2 T)- 1 -^ (2.3.1) 

+ h^A^T - 0 2 + /3/! 1 / 2 Cr 2 (l + h 1 ' 2 Ap, T ) ( - (Cr - 0 + \^ 2 Cr 2 ) , 

in L 2 ^(0, h~ s ), a dr^j . Here Ap iT is a function of (/ 0,t ) and satisfies, for all h E (0,1), 

\Ap, r \<M(p + l)T. 

The domain of the operator p ^h is 

Llfn^p^ ={ueH 2 ({0,h~ 6 )) : u\0) = -u(0) and u{h~ 5 ) = 0}. (2.3.2) 

The operator H^ t p^ : h is the Friedrichs extension in L 2 (^(0,h~ s );adT^j associated with the qua¬ 
dratic form 


qc,p,h,d u ) = 



(l + h 1/2 Ap jT ) (V(l 


\ph l / 2 T) - 2 ) (1 -(0+mh a )h 1 / 2 r) dr-|w(0)| 5 


The operator 'H^pgh. i s with compact resolvent. The strictly increasing sequence of the eigen¬ 
values of / H(,p,^,h is denoted by (A n {'H<:,p£,h))nm- 


2.3.1. Harmonic oscillator on an interval. Here we study the operator in (12.3. ip for 0 = 0, m = 0 
and Ap r = 0 which becomes the harmonic oscillator 

=-^2 + (Ct-0 2 in L 2 ((0,h~ 5 y,dr) , (2.3.3) 

and with the boundary conditions u'( 0 ) = —it( 0 ) and u(h~ s ) = 0 . 

By comparison of the quadratic forms of the operators H^p^h and %£,(),£,ft) we get that the 
spectrum of p % t h is localized near that of p £ t h as h goes to 0. This gives us a rough 
information about the spectrum of the operator p^h precisely stated in: 

Lemma 2.7. Let 0 < ci < c^, 0 < e < | and 6 E (0, 1). There exist two constants C > 0 and 
ho E (0, 1) such that, for all 

h E (0, ho), c\h e < ( < C 2 h e and \/3\ + m < C 2 , 


it holds the following. 

(1) A 2(U CM , h ) > -C\0\h^ 5 . 

(2) If e = j, 5 < \ and |£| > (3c2 + 2)h^~ 5 , then 

Ai {'H'C,p,£,h) > —1 + h2~ 2S . 

(3) Let Ao be the universal constant in Lemma \2.S[ If e < ^, 6 < ^ — 2e and |£| > AoC, then 

Ai(%£ ,p,£,h) > -i + C 2 • 

Proof. 

Step 1. 

There exists a constant C > 0 such that, for all u E i? 1 ((0, h~ 5 )), 


QCT,h,d u ) 


qc,h=o.h,d u ) ^ C\p\h2 


-<5 


qo,h(u) +IM 


L 2 (( 0 ,/i-«);dr)) 


’ U llL 2 ((0,/i-' 5 );(l-/3/i 1 /2 T ) ( i r ) 


m 


2 

L 2 (( 0 ,h~ s );dT) 


< \p\h? 


s 


L 2 ((0,h~ s );dT) ■ 
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The min-max principle yields, for all n € N, 


~ A n C%c.o,£,fc-) < C\/3\h* 5 QA rl (7^^ j o,^,fe)| + 1 


(2.3.4) 


Since the form domain of the operator 'Hharm(C)0 contains that of the operator 'Hc.fi £,h (cf. 
( 12 . 2 . 11 ) 1 . then the min-max principle yields 


AnC^c.o,?,ft) > A n (Hh arm (C.O) ■ 


(2.3.5) 


In particular, for n = 2, Lemma [2.II gives us the statement in the first item of Lemma [2.71 
Step 2. 

We estimate the quadratic form for the operator Hfi 0 , 5 ,ft as follows, for all u € H 1 ( 0, h~ s ), 

r h ~ s / 1 N 

q(,ofih{u)> J [\u , {t)\ 2 + (-2c 2 h^~ s ^ + ^ 2 )\u\ 2 JdT-\u(0)\ 2 . 

The min-max principle and Lemma 12. II yield. 

Ai(%c,o,£,h) > -1 + (|£| - 2c 2 hi~ s )\£\ . 

We insert this into (12.3.41) . That way, for |£| > ( 3 c 2 + 2)h*~ S , we get the conclusion in the second 
item of Lemma [2~n 
Step 3. 

Using (I2.3.4[) and (I2.3.5|) for n = 1, we get, 


Ai (H^fififi) > Ai (H^ t o^fi) — Ch 2 5 . 


Now, we assume that |£| > vloC) e < \ and <5 < ^ — 2e. By applying Lemma 12.21 we get, for h 
sufficiently small, the statement in the third item in Lemma 12.71 □ 


2.3.2. Lower bound for the principal eigenvalue of the operator H^fififi. In the next two propo¬ 
sitions, we determine refined lower bounds of the eigenvalue Ai (Hrfififi)- The bound is valid as 
h 0 + and is uniform with respect to the parameters £, ( and (3. 

Proposition 2.8. Let 0 < ci < C 2 , <x € (0,1), e €E (0, \), 5 € (0, \ — 2e) and n € N be the 
smallest positive integer such that 

(2n + 2)e > - . 

There exist constants C > 0 and h 0 G (0,1) such that, for all 

h € (0, ho) , ah e <(< c 2 h e , £eM, \/3\h 5 < - , \/3\ + m < c 2 , 

O 

it holds, 

Ai(« C fi£,h) > “I + e„(C) - Ph 112 - Ch r , 

where 

r = min (( 2 n + 2 )e, - + 2e, - + a) . 

Proof. Let Aq be the universal constant in Lemma [2.21 In light of the results in Remark 12.61 and 
Lemma EH the lower bound in Lemma 12.81 holds true for |£| > Aq^. It remains to prove the 
lower bound for |£| < ^4oC- 
Consider the function 

f{r) = x( r h S ) w n (r ), 

where w n {r ) is the function in Theorem 12.31 and x € C'^°([0,oo)) satisfies 

0 < y < 1 in [ 0 , 00 ), x = 1 in [ 0 , 1 / 2 ) and x = 0 in [ 1 / 2 , 00 ). 

Clearly, the function / is in the domain of the operator H^fi^fi. It is easy to check that, 


< c c 2 . 
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In light of Theorem 12.31 and the expression of p £ hf in 112. 3. ID . we may write, 

IK^C \P£,h ~ (“I + /n(C)0 “ Ph 1 ^ 2 )}f\\L 2 ((0,h- p );(l-l3h 1 / 2 T)dT) - C(( 2n+2 + h^ +2e + h^ +a ) 

< Ch r . 

(2.3.6) 

By the spectral theorem, we deduce that there exists an eigenvalue A p £,h) °f W(;,p,£,h such 
that 

WcM,h)-(~l+fn((,0)\ <Ch r . 

Now, Lemma m tells us that 

Ai ('Ht,p4,h) = ,p,t,h) ■ 

Finally, by dehnition of e n (£) in (|2.2.5[) . we have ,/ n (C, £) > e n (£). C 

Proposition 2.9. Let 0 < ci < C 2 , er € (0,1) and 5 € (0, |). There exist constants C > 0 and 
h -0 € ( 0 , 1 ) such that, for all 

1 


it holds, 


where 


h € (0, ho) , ci/i4 < £ < c 2 /i4 , (ei, |/3|h 5 <-, |/3| + m < c 2 

O 

Ai(ft c ,<8,*,ft) > - 1 + jC 2 - ^ 1/2 - , 


r = min (^1 — 46, — + crj . 

Proof. The lower bound in Lemma 12.91 trivially holds when |£| > ( 3 c 2 + 2 )hA~ s thanks to 
Lemma o 

Now we handle the case where |£| < ( 3 c 2 + 2)h*~ s . Let w\ be the function constructed in 
Theorem 12.31 and choose x £ C^QO, oo)) such that 

0 < % <1 in [ 0 , oo), x = 1 in [ 0 , 1 / 2 ) and x = 0 in [ 1 / 2 , oo). 

Consider the function 

f(r) = x(rh s )w n (r). 

Clearly, the function / is in the domain of the operator H^p^h an d 


li 2 ((0,/i-P);(l-/3/i 1 / 2 T)d r ) 


< C( 2 . 


Inserting the estimates in Theorem 12.31 into the expression of p,^,hf in (12.3.ID . and using that 
( = 0{h l / A ) and £ = 0(h*~ S ), we may write, 

IK^C ,P£,h - (—1 + /l(C)0 “ /i^ 1 ^ 2 )} / II L 2 ((0,/i - ^);(1—/3/i 1 / 2 r)dr) < C (C 4 + ^ + (C 2 + £ 2 + h a )h^ + h) 

< Ch r . 

(2.3.7) 

Now, the spectral theorem and Lemma 12.71 yield 

= -1 + /i(C,0 - Ph 1/2 + C?(/i r ). 

Noticing that min /l (C; 0 = t£ 2 > we finish the proof of Lemma [2.91 □ 

kl<(3c 2 +2)^-^ 4 
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3. Analysis of the semi-classical Laplacian with a weak magnetic field 

3.1. Semi-classical Laplacian with weak magnetic field. We will introduce a new semi- 
classical magnetic Laplacian but with a Robin condition not involving the parameters a and 7. 
These two parameters will be absorbed by a new (small) parameter £. 

For h > 0 and ( > 0, we introduce the operator 

Ph ,c = -(W-<A 0 ) 2 in L 2 (Q), (3.1.1) 

whose domain is 


D(P hX ) = {u € Hl- HAo ) : -(W-iCA 0 ) 2 €t 2 (n) and 1 /■ (fcV - iA 0 )ti 
This operator is defined via the quadratic form 


—h}/ 2 u on <9fl}. 

(3.1.2) 


u ^ q h ,d u ) = [ |(hV — i£Ao)ri | 2 dx — h 3 ^ 2 f \u\ 2 ds(x). (3.1.3) 

Jn Jon 

Let a(Ph t () be the spectrum of the operator Ph,(- We introduce the ground state energy, 

Ai(h,C) = inf cr(P hjC ). (3.1.4) 

There is a relationship between the ground state energies in (11.1.71) and (13.1.41) displayed as 
follows: 

7 4 / h 2 ~ 2a h l ~ 2a \ 

= (3.1.5) 

Now Theorem 11.11 follows from: 

Theorem 3.1. Let 0 < c\ < C 2 and e > 0. Suppose that 

0 < h < 1 and c\h e < £ < C 2 h e . 


It holds the following. 

(1) If e < then as h —>• 0+, the ground state energy in (13. 1 .4[) satisfies 

Ai(h,C) = -h + e n (C)h - «: max /i 3/2 + h 3 / 2 o(l), 

where e n (£) is introduced in (12.2.51) and n £ N is the smallest positive integer such that 
(2n T 2)e > ^. 

(2) If e = j, then as h 0 +; the ground state energy in (13.1.41) satisfies 

Ai (h, C) = ~h + ^( 2 h - K mSLX h 3/2 + h 3/2 o( 1 ). 

(3) If e > j, then as h —> 0 +; the ground state energy in (13.1.41) satisfies 

Ai (h, C ) = -h- K max h 3/2 + h 3/2 o( 1 ). 

(4) There exist constants p £ (0, ^), p* £ (0, -j), C > 0 and ho £ (0,1) such that, for all 
h € (0, ho), every ground state Uh of X\(h, £) satisfies, 

IKclU 2 (n bnd ) < Cl exp , ||^,cllL 2 ( n int) < exp 

where 

flint = {x £ fl : dist(x, <912) > h p } and 

flbnd — {% £ fl \ flint • ^max ^(p(®)) — h ^ } . 



The proof of the items (l)-(3) in Theorem 13. II follows from Proposition 13.21 and Proposition [3761 
We will give explicit bounds to the remainder h 3 / 2 o(l) in the form 0{h r ) where r depends on e 
and satisfies r > More specifically, we find that 

r = min(r*, r *), 


where r* and r* are introduced in (|3.3.2I) and (13.5.31) respectively. 
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The proof of the item (4) in Theorem 13.11 follows from Theorems 13.31 and 13.91 


3.2. Boundary coordinates. We will perform various computations of trial functions sup¬ 
ported in a tubular neighborhood of the boundary. To single out the influence of the boundary 
curvature, we need a special coordinate system displaying the arc-length along the boundary and 
the normal distance to the boundary. We will refer to such coordinates as boundary coordinates. 
These are the same coordinates used in the semi-classical analysis of the magnetic Laplcian (cf. 

EE]). 

The boundary coordinates are valid in every connected component of the boundary. For 
simplicity, we will suppose that <914 has one connected component; if more than one connected 
component exists, then we use the coordinates in each connected component independently. Let 

R/(|<914|Z) 3m M(s ) £ <914 

be the arc-length parametrization of <914 and oriented counterclockwise. At the point M(s ) £ <914, 
n(s) is the unit outward normal vector; the unit tangent vector T(s) and the curvature k(s) are 
defined as follows 

T(s) := M'(s), 

and 

T'(s) = k(s) v(s). 

The counterclockwise orientation of the parametrization is displayed as follows, 

Vs € R/(|<914|Z), det(T(s),i/(s)) = 1. 

The smoothness of the boundary yields the existence of a constant to > 0 such that, upon 
defining 

Vt 0 = {x £ 14 : dist(x, <914) < to}, 

the map 

<h : M/(|<914|Z) x (0,to) 9 (s,t) m x = M(s) — tv{s ) £ Vt 0 . 
becomes a diffeomorphism. Let us note that, for x £ Vt 0 , one can write 

$ - 1 (x) := (s(x),t(x)) £ R/(|<914|Z) x (0,t 0 ), (3.2.1) 

where t(x) = dist(x, 4414) and s(x) £ R/(|<914|Z) is (uniquely) defined via the reation dist(x, <914) = 
\x — M(s(x))\. 

Now we express various integrals in the new coordinates (s,t). First, note that the Jacobian 
determinant of the transformation 4 > _1 is given by: 


a(s, t) = 1 — t«(s). 


In the new coordinates, the components of the vector field Ao are given as follows, 

.— fj'Y 

Ai(s, t) = A 0 • — = (1 - tAt(s))A 0 (4>(s, t)) • M'(s), 

(IT 

A 2 (s,t) = A 0 • — = A 0 (<F(s,t)) • v(s). 

The new magnetic potential Ao = (Ai, A 2 ) satisfies, 
dA dA i 

2 (s,t) -— ~(s,t) ds A dt = curl Ao(4> - 1 (s, t))dx A dy = (1 — tn(s))ds A dt. 


(3.2.2) 


ds dt 

For all u £ L 2 (Vs), we assign the function u defined in the new coordinates as follows 

u(s, t) := u(&(s, t)). (3.2.3) 

Consequently, for all u £ H l (Vt 0 ), we have, with u = u o$, 

|(W — i(Ao)u\ 2 dx = J (1 — tK,(s))~ 2 \(hd s — i(Ai)u\ 2 + | (hdt — iCA 2 )u| 2 (1 — tn(s))dsdt , 

(3.2.4) 
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and 


/ \u\ 2 dx = / \u(s,t)\ 2 (l — tK,(s))dsdt. 

Jv t() J 

/ \u\ 2 dx = / |u(s, f = 0)| 2 ds. 

dv tn nan J 


(3.2.5) 


(3.2.6) 


Finally, we recall a useful gauge transformation that we borrow from um- Let xo €E <912 and 
V xo be a neighborhood of xo in 12. There exists a smooth function <f XQ in <h^ 1 (V X0 ) such that, in 
the boundary coordinates, 

, +2 

(3.2.7) 


A ~ V( s ,t)4>xo = ( - t + itk(s), 0 


3.3. Upper bound for the principal eigenvalue. In the rest of this paper, we will use the 
following notation. For all e > 0 and For all f € (0,1), define 


MC) = 


e n (C) if e < 1/4, 
K 2 if e = 1/4, 
0 if e > 1/4, 


(3.3.1) 


where n € N is the smallest positive integer satisfying ( 2 n + 2 )e > 77 . 

Proposition 3.2. Under the assumptions in Theorem \S. 1[ there exist two constants C > 0 and 
ho € (0,1) such that, for all h € (0, ho), the ground state energy in 113.1.411 satisfies, 

Ai (h, ()<-h + bfi()h - Kmax h 3 / 2 + Ch r *, 

where 

_j 1 + min ((2n + 2)e, \ + e) if e < 1/4, ( 332 ) 


r* = 


13/8 


if e > 1/4. 


Here b e (() is as in (13.3.111 . 


Proof. The proof consists of constructing a trial function Vh(x ) and computing its energy. This 
trial function will be defined via the boundary coordinates (s,f) in (13.2.111 . Select xq € <912 such 
that 

k(s(x 0 )) = K max 

is equal to the maximal curvature. We may choose the coordinates ( s , t ) in (13.2. ID such that 
s(x 0 ) = 0. Let V xo be a neighborhood of the point xo in 12, and tpo = (p XQ be the function defined 
in V XQ and satisfying (I3.2.7D . 

The construction of the trial function Vh and the computation of its energy will be done for 
the cases e < e = | and e > \ independently. 

The case e < 

Let n € N be the largest positive integer such that (2 n + 2)e > /. Recall the definition of 
fn( C,0 i n (|2.2.6D . Select = £ n (C) such that 

fn(C,€n) = e n {C) = min{/ n (C,0 : |£| < Cmax(A 0 ,1)} and |£ n | < Cmax(A 0 ,1). 

The trial function Vh is defined using the (s, ^-coordinates and the relation in (13. 2. 3D as follows, 


l+e 


v h {s,t) =ch 4 Xi 


Xi 


fc/ 2 


~ lipo w n (h 1/2 t ) exp 


(3.3.3) 


Several objects appear in the definition of vy t : 

(1) xi € C/°(M) satisfies 0 < x < 1 hi K, xi = 1 hi [—1/2,1/2] and suppxi C [—1,1] ; 

( 2 ) c = lixi llx ,2 (K.) ’ 

(3) w n {r) in the function in Theorem 12.31 : 

( 4 ) P= 3 + e . 
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The upper bound in Proposition 13.21 follows from the min-max principle and the following two 
estimates: 


\\v h \\mn)-l\ <Ch x l\ 

{Ph,C + h — e n {Q + Kmax^ 3 / 2 )lhj||L 2 (Q) < h r * , 


(3.3.4) 

(3.3.5) 


where r* = 1 + min ^(2n + 2)e, \ + is given in (13. 3. 2D . 

The estimate in (I3.3.4P is easy to obtain in light of the expression of Vh and the formula 
in (I3.2.5p . For the estimate in (I3.3.5D . notice that, after expressing the operator P^.h in the 
boundary coordinates (s,t), we get (compare with (I3.2.4p l 


{Ph, C T h Gn{C) T Kmaxfr ^ ||z, 2 (f2) — + h Sn(C ) T K max /l ^)uh\\l t '2^ a dsdt) j 


where 


a(s, t) = 1 — tn(s) 


1+e f t 


u h (s,t) =ch 4 xi 


,{h 1/2 t ), 


and 


Note that, 


where 


and 


hP J A1 W/ 2 - 

L h ,c = ~h 2 cT 1 d t (ad t ) - a~ 2 [hd s - i(t(l - ^k(s) 

^h,c, u h = chr^xi PL hX xi w n( hl/2 t ) + R h, 

PL h ,c = -h 2 aT 1 dt(ad t ) + a~ 2 (c* (l - ^«(s)) - /i 1/2 £n) , 


(3.3.6) 

(3.3.7) 


R h = ch 1+ * e Xl f1- ) Wn (h 1/2 t) h 2 d 2 xi 


, -i+i f t 
+ 2icb 4 Xl l — 


w. 




h e / 2 Jl 

(h 1/2 Zn-(t(l 


2^ S 


hd< 


Xi 


h*! 2 


(3.3.8) 

(3.3.9) 


It is easy to check that 

11 R h 11 L 2 (adsdt) < + Ch 2 " e < C7li+f . 

Now we perform the change of variable t = h l ^ 2 T and get (a = 1 — /i 1//2 tk(s)), 

r / h 1 / 2 r \ 2l 

PL hX = h - a~ 1 d T (ad T ) + a~ 2 ^t( 1-re(s)) - £ n j 

where the operator 'H^ tK ( s )Xn,h i s introduced in (12.3. ID with /3 = k(s), A, g jT = /i -1 / 2 (a -2 — 1), 
m = 0 and <5 = -7 — p = -(77 — 2e) < 77 — 2e. Now, it is easy to prove that (compare with (j2.3.6p ). 


(h 1 PL hX - ( - 1 + / n (C,Cn) - «(s)h 1/2 )^Xi (h2 p T^W n (T) 


L 2 (adr) 


< Ch r , 


where 


r = min ((2 n + 2)e, - + 2e) . 


Returning back to the t variable then integrating with respect to the s variable, we get, 


{PLfi, C~ (- h + h 1/2 fn(C,£n) - K(s)h 1/2 )^X1 (h p t)w n (h 1/2 t) 


L 2 (adsdr) 


< Ch* +r . 
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Now, we insert this and (13. 3. 91) into (13. 3. 61) to get, 


{L h ,c + h-e n (0 + K(s)h 3 / 2 )u h \\ L 2 iadsdt) <Ch = Ch r * , 

To finish the proof, we notice that K max = /c(0) and in the support of the function Uh, we have 

|k(s) -«(0)| < Ch 3 ' 8 . 

The case e = 

Now the trial function Vh is defined using the (s, ^-coordinates and the relation in (13.2.31) as 
follows, 


v h (s,t) = ch 5/16 Xi (j^j Xi 


h 1 / 8 


exp (* 


(3.3.10) 


where the constant c and the function xi are as i n (13.3.31) . and w\ is the function defined in 
Theorem 12.31 

Performing a calculation similar to the one done for the case e < 4 (in particular, using (I2.3.7P 
for 5 = jg and f3 = /c(s)) we get, 


\\ v h\\L*(n) ~ 1 


< Ch 1/2 , 


{Ph£ + h C ) 2 + Kmaxh 8 ^) 


Vh 




< Ch 13 / 8 


The min-max principle now yields the desired upper bound for X\(h, £). 

The case e > 

Here we simply take the same trial state for the case without a magnetic field but times a 
phase (cf. [7] [12]). Precisely, we define Vh as follows, 


Vh 


( s,t) = ch Xi (^g) e !W !*o(ft 1/2 t), (3.3.11) 

where the constant c and the function xi are as hi (I3.3.3P and uq(t ) = \/2 exp(—r). Easy 
calculations similar to those done for e < -j give us (cf. [7. 12]) 

IKIli^n)-! <Ch 1/2 , 

{Ph,c + h + K mSLX h 3/2 )v h \\ L 2^ < h 13/8 . 

The min-max principle now yields the upper bound for Ai (h,£). □ 


3.4. Concentration of bound states near the boundary. 

Theorem 3.3. Let 0 < ci < c^, e > 0 and a < 1. There exist constants C > 0 and ho £ (0,1) 
such that, if h £ (0, ho), C € (c\h e ,C 2 h e ) and uj,^ is a L 2 -normalized ground state of Ph.Q, then, 

J ^ (\uh,s(x) I 2 + h _1 |(/iV - iCAoWxOr)! 2 ) exp dx<C . 

The proof of Theorem 13.31 makes use of the result in: 

Lemma 3.4. Under the assumptions in Theorem \3.tA if 0 < p < w £ iP 1 (H) and supprc C 
{dist(x, <9fl) < 2 h p }, then 

f \(h\7 — i(Ao)w\ 2 dx — h 3 P f \w\ 2 dx > — ^ f \w\ 2 dx . 

Jn Jon 2 Jn 
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Proof. Let qh,((') be the quadratic form in (13.1.31) . The diamagnetic inequality yields, 

qh,(( w ) — [ | hV|re| | 2 dx — h 3 / 2 f \w\ 2 ds(x). 

Jn Jan 

In boundary coordinates, the inequality reads (cf. (|3.2.4l) h 


Qh,c( w ) — (1 — ChP) JJ | hVv | 2 dsdt — h 3 ^ 2 j |u( 


s. t = 0)| 2 ds . 


where v = |u;o3>(s,i)| (cf. (I3.2.3p l. Applying the change of the variable t = hJ^T and comparing 
with the operator in (12.1. ID . we get, 


□ 


qh,((w) >—(1 — Ch p )h\iCHofi) // |u(s, t)\dsdt = — (1 — Ch p )h // \v(s,t)\dsdt. 


Returning back to Cartesian coordinates, we get the inequality in Lemma 13.41 

Proof of Theorem \3.JA The proof is similar to that of Theorem 5.1 in [7J. Let t(x) = dist(x,cAl) 
and <L(x) = exp( °*^ ). We perform an integration by parts to write the following identity, 

qh(uh X )-= [ {\(hV-i(A 0 )($u hX )\ 2 -h 2 \V$\ 2 \u hX \ 2 ) dx-h 3 ' 2 [ \<S>u h , c \ 2 ds(x) 

Jn Jan 

= \i(h, C)||$ ^/i,cIIl2(q) • 

Consider a partition of unity of R 

xi + xl = i, 

such that xi = 1 hi (—oo, 1), suppxi C (—oo, 2), xi > 0 and X 2 > 0 in R. 

Define 


(3.4.1) 


Xj,h ( x ) = Xj » j e {!. 2 } 


Associated with this partition of unity, we have the simple standard decomposition 

2 

?h,cKc) = J2 q thx( Uh x) > 

3 = 1 

where (by Lemma. l3.4p 

9gh,c( u Ac) = [ {\(hV-i(A 0 )(xi,h$Uh,c)\ 2 - h 2 \V(xi,h$)\ 2 \uh,c\ 2 ) dx 


(3.4.2) 


~h 3/2 f \xi,h$ ^fe.cl 2 ds(x) 

Jan 


(3.4.3) 


>-77 [ \xi,h$u hX \ 2 dx-Ch [ \<$> u h x\ 2 dx >-Ch . 

1 Jn Jn 


and 


(3.4.4) 


Q 2 ,h,d u h,d= / (l(hV-iCA 0 )(x 2 ,i t d > w/ l ,c)| 2 -h 2 |V(x 2 ,/ l ^)| 2 |ibi,cl 2 ) dx ■ 

Jn 

The definition of <h and the fact |Vt(x)| = 1 a.e. together yield 

/ {X 2 ,h$)\ 2 \uh ,(\ 2 dx < a 2 h / \x 2 ,h$u h ,<;\ 2 dx + Ch / |ri^| 2 cte. 

J q J J 

We insert this and (13.4.31) into (I3.4.ip . write \\{h, Cf) < |(—1 — a 2 ) by Proposition 13.21 and 
rearrange the terms to obtain, 

(j(hV - i(Ao)(x 2 ,h & u hx)\ 2 + - n 2 )h\X 2 ,h $ Uh^l 2 ^ dx < Ch. (3.4.5) 

This is enough to deduce the estimate in Theorem 13.31 □ 
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We record the following simple corollary of Theorem 13.31 

Corollary 3.5. Let (0,|), e > 0 and 0 < ci < C2- There exists ho E (0,1) such that, for 
all h E (0, ho) and ( € (c\h e ,C2h e ), every L 2 -normalized ground state u^ of the operator 
satisfies, 

/ \uh,(\ 2 dx < ex.p (—cih p ~2) . 

J cihP<dist(x,dfl)<C2hP ^ ' 

3.5. Lower bound for the principal eigenvalue. 


Proposition 3.6. Let e > 0 and 0 < ci < C 2 - There exist constants C > 0, ho E (0,1) and 
r* > | such that, for all h E (0,/io) an d C £ {cih e ,C 2 h e ), the ground state energy in (13.1.41) 
satisfies, 

Ai (h, ()>-h + b e {()h - K max h 3 / 2 - Ch r *, 
where 6 e (C) is introduced in (13.3.11) . 

For e > i Proposition 13.61 follows from: 

Lemma 3.7. Suppose that e > j. Under the assumption in Proposition [X5] for all u in the 
form domain of the operator Ph£, 


Qh,c(u) > / U h ^{x)\u\ 2 dx . 


where 


TT , . f — h — n(s(x))h 3 / 2 — ChJ/ A if disthr, dUt) < 2/j 1 / 8 , 
t, « W = |o tfdtot(z,8S2)>2h'/ 8 , 

and qh,((-) is the quadratic form in (13.1.31) . 

Proof. This is a consequence of the diamagnetic inequality and [7J Thm. 5.2]. 


□ 


In the case e < Proposition 13.61 is a consequence of Lemma [3.81 below (applied with w = u ^ 
and Uh.( a L 2 normalized ground state of the operator Ph.cf) and the variational min-max principle. 
The constant r* in Proposition 13.61 depends on e. It is introduced as follows. For all e > 0, let 

( l min (2e, 1 - 4e) if e < j , 

1 1/8 if e > 1, 

| — \ min (2e, 1 — 4e) if e < 1/4 
7/16 if e = 1/4, 

1/8 if e > 1/4, 

min ^1 T (fin T 2)e, ^ + 2e, ^ + a, 2e T 4 p T 2cr — 1 4 u T p, 2 — 2 cj^ 

7 
4 

where n E N is the smallest positive integer satisfying (2 n + 2)e > /. 

Note that 0<p<|,0<cr<l and when e < the following three conditions are satisfied 

2e + 4p + 2a — ^ ^ , 

2 - 2a > |, (3.5.4) 

p + a > \ . 

Consequently, for all e > 0, the number r* satisfies 

„ 3 




(3.5.1) 


(3.5.2) 

if e < 3 , 
if e > 3 , 

(3.5.3) 
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Lemma 3.8. Let M > 0, 0 < e < 4 and 0 < c\ < c^. There exist two constants C > 0 and 
ho G (0,1) such that, for all h G (0, ho) and ( G ( c\h e , c^hL), if w is a L 2 normalized function in 
the form domain of the operator P^^ and 


exp 


dist(x,<9fl)\ 

2 W ) 


w 


L 2 (n) 


<M, 


for some p < p' < b, then it holds the following: 


Here 


U h> c (x) = 


Qh,c{w)> [ U h ^(x)\w\ 2 dx — Ch 2 
J n 


_ J —h + b e ( C) — k(s(x))/i 3//2 — Ch r * if dist(x, dfl) < 2 h p , 


(3.5.5) 


(3.5.6) 


0 if dist(x, dfl) > 2 h p ; 

• a, p and r* are introduced in (13.5.1(1 . (13.5.2(1 and (13.5.311 respectively; 

• b e {(f) is introduced in (13.3.1(1 ; 

• Qh,((') is the quadratic form introduced in (13.1.31) . 

Proof of Lemma \3.tA 

The lengthy proof of Lemma [3.81 is divided into four steps. 

Step 1. Localization near the boundary. 

Consider a partition of unity of M, 

xi + x! = i 

with Xi = 1 in (— oo, 1], suppyi C (— oo,2] and supp %2 C [l,oo). For j G {1,2}, put, 


Xj,h{x) = Xj ( 


( dist(x, dfl) 


hP 


We have the decomposition 

Qh,d w ) = qh,c(Xi,hw) + Qh,g{X2,hw) - h 2 ^ || \Vxj,h\w ||^ 2(n) 

where 


l=i 


and by (13.5.511 . 


Thus, 


Qh,dX 2 ,hw) = / \hX7( X 2,hw)\ 2 dx>0, 

■In 

h 2 1| \^Xj,h\w II Ja (n) < h 2 ~ 2p exp = 0(h °°). 


(3.5.7) 


Qh,d w ) ^ Qh,dxi,h w) + 0(h°°). 

Step 2. Analysis near the boundary. 

Let us cover the boundary dfl by a family of open disks ( B(xj , /i 1 / 8 )). Let ( fj) C C°°(dfl) be 
a partition of unity in dfl such that, for all j, 

supp fj C B(xj, h a ) n dfl , and \V fj\ < Ch~ a . 

We extend fj in the tubular neighborhood {dist(x, dfl) < 2 h p } of the boundary via the formula 

fj( x ) = fj( s ( x ))- 

We decompose the boundary term in (13.5.711 as follows, 

Qh,dXl,hw) = ^q h ,dfjXl,hw) ~ h 2 ^ II |V/j|xi,ftw||ia (n) 


> 


5 Zflh,dfjXi,hW) - Ch 2 2a \\xi,hw\\h(n) ■ 


(3.5.8) 
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We will write a lower bound for each term qh,((fjXi,h w ) as follows. First, let us denote by 

Kj = K(s(xj)) . 

By smoothness of the scalar curvature and boundedness of the boundary, we know that 

|k(s(x)) — Kj | < mh a in B(xj, h a ), 


where 


m = sup |k , (s(x))| . 
x£dQ 


That way we get the following pointwise lower bound in every B{xj , h a ), 


| (hd s — i(t( 1 — Jk(s))w \ 2 > (1 — h l ^ 2 )\{hd s — i(t(l — JtKj)w \ 2 — 4i( 2 i 4 /i 2<T 2 \w \ 2 . 


Let (f)j = (f> Xj be the function satisfying (13.2.711 in B(xj,2h a + 2 h p ). Dehne the function 
Dj = fjXi.hW . We express the quadratic form qhx(fjXi,h w ) l n boundary coordinates and 
then we use the aforementioned inequalities to write, 


qh,({fjXi,h,w) > 

JJ (\hdtVj\ 2 + (1 — h 1 ^ 2 )(l — tKj — Ch a t)~ 2 \(hd s — i(t(l — ^tKj)w\ 2 ^J (1 


Kjt — mh a t)dsdt 


J \vj(s,t = 0)\ 2 ds — C 


( 2 t 4 h 2a 2 |uj| 2 (l — Kjt — mh a t)dsdt. 


Let 


(3.5.9) 


5 = ^ — p, f3 = Kj, A p T = h 1//2 (1 — h 1//2 )(l — (Kj + C7i <t )/i 1//2 t) 2 —1 


(3.5.10) 


Note that, for t € (0, h^~ p ) and r = h~^t, |A^ T | < C(\f3\ + 1 )r. Thus, we can apply the results 
in Sec. 12.31 

We return back to (13.5.911 . Note that Q = 0(h e ) and in the support of Vj , the term t 4 is of 
order 0{h Ap ). We apply the change of variable t = h}^ 2 r then the Fourier transform with respect 
to the variable s to obtain, 


qh,dfjXi,h,w) > JJ {(fcgAi^c^))-^ 2 ^ 2 1 \f jX i,hw\ 2 (l - Kjt - mh a t)dsdt. 

(3.5.11) 

Step 3. Lower bound in the case e < 

By the assumption on p and a, we find that 

0 < <5 =- p < -2e 

2^2 

so that we can apply Proposition 12.81 Let r* > | be the constant introduced in (I3.5.3P We infer 
from (13.5.1111 . 


qh,dfjXi,hw) > JJ |-h + /te n (C) - Kjh l/2 -Ch r *}\fjXi,hw\ 2 (l -tKj -mh a t)dsdt. (3.5.12) 

Now, in (13. 5. IIP , we replace Kj by k(s) + 0(h a ) and use that a + p > -7 to replace the term 
1 — tKj — mh a t by 1 — tn(s) and get, 

qh,dfjXi,hw)> JJ j-h + he n (() - K{s)h 1/2 - C7t r *} \fjXxhw\ 2 iX ~ tn(s))dsdt . (3.5.13) 
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We insert (13.5.1311 into (13. 5. 811 and obtain 


Qh,c(Xi,hw) > ^2 



-h + he n ( C) — K{s)h 1 ^ 2 — Ch r * j |/jXi,^ui| 2 (1 — tn(s)) dsdt — Ch 2 

= j j —h + he n ( C) — k(s(x))/i 1//2 — Ch r * | |xi ,h w \ 2 dx — Ch 2 . 

Using that \xi,ti w \ < |io| and that — h + he n (() — k(s)/i 1//2 — Ch r * < 0, we get further, 

Qh,dxi,h w ) = J |-h + he n (C) ~ n(s)h 1/2 - Ch r * j \w\ 2 dx - Ch 2 . 

Finally, we insert this into (13. 5. 71) to get (13.5.6[) . 

Step 4. Lower bound in the case e = 

The analysis here is similar to that in Step 3 and we will be rather succinct. Note that the 
assumption on 5 and p ensure that 0 < <5 = ^ — p = ^ < | so that we can apply Proposition 12.91 
and infer from (13.5.111) . 


qh,dfjXi,hw) > J J j-h + j( 2 h - Kjh 3/2 - Ch r *| \fjXi,h'w\ 2 (l - Kjt - mh a t) dsdt 

h + ~ k(s)/i 3 / 2 — Ch r | |/jXi,/iH 2 (l — K{s)t) dsdt. 

We insert this into (13.5.811 and (13. 5. 71) to get ([3.5.611 for e = 


(3.5.14) 


□ 


3.6. Concentration of ground states near the points of maximal curvature. 

Theorem 3.9. Let 0 < ci < C2 and e > 0. There exist constants p G (0, ^), p* € (0, -j), C > 0 

and ho G (0,1) such that, for all h G (0, ho), C € (cih e , C 2 h e ) and u^ a normalized ground of the 

operator Ph.q, 

f \ u h,(\ 2 dx < Cexp (—h v *^2 . 

J{dist(x,dQ)<h q }r\{Kmax.—K{s(x))>h r i*} ' ' 

We will prove Theorem 13.91 in the case 0 < e < The case e > ^ is a standard consequence 
of the inequality in Lemma 13.71 (cf. [2] Thm. 8.3.4]). 

An important ingredient in the proof of Theorem 13.91 is: 

Lemma 3.10. Let e > 0 and p be as in (I3.5.2[) . Let p' G (p, ^). There exist two constants C > 0 

and ho G (0,1) such that, for all h G (0, ho) and u in the form domain of the operator Ph£, 

Qh,<:( u )> / Vh,d x )\ u \ 2 dx , 

Jn 


where 


, . _ f —h/2 if dist(x,<9fl) > 2 h p , 

h ’ C ‘ X \ —h + b e (()h — K,(s(x))h 3//2 — Ch r if dist(x, d£l) < 2 h p , 


and qh.Q is the quadratic form in (13.1.31) . 

Proof. Let ft be the ground state energy of the operator Ph.g — Vh^- We will prove that ft > 0. 
The min-max principle and Theorem 13.31 together yield 

ft < ((Ph, c - Vh,c) u h,c, u h,(;)L 2 (n) = W(h, C) - / Vh,g\uh,c\ 2 dx < Ch 7 . 

Jn 

Let w he a L 2 normalized ground state of the operator P^x ~ Vh,C- We will prove that, 

t(x 


exp 


hP' 


l 2 ( a) 


<C, 


(3.6.1) 
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where t(x) = dist(x, <9fi) and 4>(x) = exp(^r). We perform an integration by parts to write the 
following identity, 




\uhx)'-= [ (|(W - i(A 0 )(<f>w)\ 2 - V h \$w\ 2 - h 2 \X7$\ 2 \u hX \ 2 ) dx - h 3/2 [ |$u;| 2 ds(x) 
Jn ’ Jan 

= V\\$W\\ 2 L 2 (n) . 

Consider a partition of unity of R 

xl + xl = 1 , 

such that xi = 1 in (—oo, 1), suppxi C (—oo, 2), xi > 0 and %2 > 0 in R. 

Define 


(3.6.2) 


Xj,h(x) = Xj (^Pj , j£{ 1,2} 


Associated with this partition of unity, we have the simple standard decomposition 

oft, cH > 


(3.6.3) 


3 = 1 


where 


= / (l(^ v - *C A o)(xi,ft,$w )| 2 - V h £\xi,h$w\ 2 ~ ^ 2 |V(xi,h^)| 2 |^| 2 ) dx 
J n 

-h 3/2 f \xi,h$ w\ 2 ds(x) , (3.6.4) 
Jan 

and 

Q2h,c( w ) = [ (K^-V-iCA 0 )(x 2 ,h^^)| 2 - V ht(: \x 2 ,h^w\ 2 - h 2 \X7(x2,h$)\' 2 \w\ 2 ) dx . (3.6.5) 

Jn 

Lemma [3.41 the bound Vh£ < 0 and the normalization of Uh £ together yield 

Qi,h,d u ^0 - ~ Ch - 

We insert this into (13.6.31) . then we insert the obtained inequality into (13.4.11) . use the bounds 
— Vh£ > 7jh, p < Ch r * and then rearrange the terms to obtain, 

^\(hV+ iCA 0 )( X 2 ,h^u hX )\ 2 + h^-Ch 1 - 2 P' -C'/i r *- 1 )|x2,/ l $«h,cl 2 ) dx <Ch. (3.6.6) 

Using the inequality ^ — Ch 1_2 p' — Ch r *~ 1 > | then dividing by h, we get 

j ^ ^/i _1 |(/iV - iCA 0 )(x 2 ,/ l ^^)| 2 + ^|X 2 ,fc$w| 2 ^ dx < C. 

This is enough to deduce the estimate in (13.6.11) . 

Having proved (13.6.11) . we may use the result in Lemma [3.81 and write, 


Qh,<;(™)> [ U h ^{x)\w\ 2 dx - Ch 2 
Jn 


Note that, by selecting C > 0 sufficiently large, we may write 

U h x(x) - Ch 2 > V h:( (x) in D . 
This finishes the proof of Lemma 13.101 


□ 
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Proof of Theorem \S. .91 Define the function i f>(s) = A max — k(s). The function cf defines a C 1 
function in dQ. The boundedness of 00 ensures the existence of a constants Co > 0 such that 

V s € (0, |00|), \<p\s)\ 2 <Ccf(s). 

Let a and p be as in (13. 5. ID and (13. 5. 2D . Choose x € C'^°(R) such that 0 < x < 1 hi M, x = 1 
in [to To] and suppx C [—1,1]. Here to < 1 is a geometric constant such that the boundary 
coordinates (s,t) are valid in the tubular neighborhood {dist(x,00) < to} (cf. Sec. 13.2D . 

Define the function 

-t(x) = exp 

where t(x) = dist(x,00) and <5 G (0,1). We will fix a choice for 5 at a later point. 

Let us write the following decomposition formula 

^i(h, C) ||^2 (q-) + h / 5 = Qh,( i^ u h,() > 

where 

${x) = x{t{x))(t>{s{x)). 

Using Theorem 13.31 and the bound \cj)'\ 2 < C(f, we may write, 

W^^u^Wl^^Oih^ + lC [[ l^smuh^dsdt. 

We use the upper bound for Ai(6, £) in Proposition [321 the lower bound for qhx(') in Lemma |3.101 
and the simple lower bound ^ — Ai(6, () > \h in { t(x ) > 26 p }, we get, 


t<2hP 


s ) - 2 Ch 3 / 2 5 2 cf(s) - C 6 r ‘) \^u hX \ 2 dsdt < CK 


We choose 6 = and get 


t<2t 0 


— Ch r 2 ] \^uhx\ 2 dsdt < Ch r . 


In particular, setting p = max(r* — §, 7 ), we write 


t<2h 2p 
4>(s)>4Ch' n 

Selecting r/* € (0, rj) finishes the proof of Theorem 13.9 


<f>(s)\'f>Uh t (\ 2 dsdt < C. 


□ 


4. Analysis of Diamagnetism 

Proof of Theorem II. 31 There exists a simple relationship between the eigenvalues in (II. LTD 
and (ll.2.2[) . This relationship is displayed as follows 

MP;H) = fL 2 /ri(h;6,a,7), 

where h = H _1 , 6=1 and 7 = j3H~ 1+a . 

The assumption in Theorem 11. 3 1 ensure that 

• As /3 —>• — 00 , the semi-classical parameter h —>• 0+ ; 

• The parameter 7 is uniformly bounded, i.e. 7 = 0(1) as (5 —>• — 00 . 

The ground state energy pi(h] b, a, 7 ) is estimated in Theorem 11.11 for a < | (more precisely, 
this is a consequence of Theorem im This yields the estimates announced for /2i(/3;fL) in 
Theorem 11.31 for a < 

For a > ^, the estimates in Theorem 11.31 follow from the following result proved in [10j| 

( 0 ( 0)6 + 6 o(l) if a > \ , 

Mi (h-,n = l,o;, 7 ) = < 

[ @( 7 ) 6 . + /io(l) if ol = \ . 
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